Thermal gradient-induced forces on geodetic reference masses for LISA 
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The low frequency sensitivity of space-borne gravitational wave observatories will depend criti- 
cally on the geodetic purity of the trajectories of orbiting test masses. Fluctuations in the temper- 
ature difference across the enclosure surrounding the free-falling test mass can produce noisy forces 
through several processes, including the radiometric effect, radiation pressure, and outgassing. We 
present here a detailed experimental investigation of thermal gradient-induced forces for the LISA 
gravitational wave mission and the LISA Pathfinder, employing high resolution torsion pendulum 
measurements of the torque on a LISA-like test mass suspended inside a prototype of the LISA grav- 
itational reference sensor that will surround the test mass in orbit. The measurement campaign, 
accompanied by numerical simulations of the radiometric and radiation pressure effects, allows a 
more accurate and representative characterization of thermal-gradient forces in the specific geometry 
and environment relevant to LISA free-fall. The pressure dependence of the measured torques allows 
clear identification of the radiometric effect, in quantitative agreement with the model developed. 
In the limit of zero gas pressure, the measurements are most likely dominated by outgassing, but at 
a low level that does not threaten the LISA sensitivity goals. 

PACS numbers: 04.80.Nn, 07.87.-l-v, 95.55.Ym 
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I. INTRODUCTION 

LISA, the Laser Interferometer Space Antenna is 
a ESA-NASA mission to create the first spaceborne in- 
terferometric observatory of gravitational waves, in the 
frequency range between 30 /iHz and 1 Hz. LISA should 
permit not only detection, but detailed, high resolution, 
long measurement time observation of gravitational wave 
signals, opening the possibility for important discovery in 
fundamental physics, astrophysics, and cosmology. 

LISA consists of three identical spacecraft separated 
by 5 X 10^ km, forming a nearly equilateral triangle that 
orbits the sun at a distance of 1 AU. Inside each space- 
craft, there are two "free- falling" cubic test-masses (TM) 
in nominally geodesic motion. Gravitational waves are 
detected by an interferometric measurement of the rela- 
tive change in the distances between distant free-falling 
TM caused by the gravitational wave metric perturba- 
tion. LISA aims at a gravitational wave strain noise floor 
of 4 X 10~^^/IIz^/^ near 3 mHz, increasing as 1/ P at 
lower frequencies. Performance at these lower frequen- 
cies is limited by the purity of the TM geodetic motion 
and requires perfect free-fall, along the sensitive intcrfer- 
ometry axis (referred as the x axis here) , to within an ac- 
celeration noise of 3 x lO^^^m/sVHzi/^ (3fm/sVHzi/2). 
For the 2 kg TM foreseen for LISA, this is equivalent to 
a force noise of 6 fN/Hz^/^. While placing test particles 
in nearly perfect geodetic orbits inside of a controlled, 
co-orbiting spacecraft is essential to various space grav- 
itational experiments, the extremely low force noise re- 
quirement and low measurement frequencies make LISA 
even more demanding in this respect than other previous 
(GPB 0) and future (STEP g. Microscope Q) mis- 
sions. 

While the spacecraft shields the TM from environmen- 



tal disturbances such as solar radiation pressure, it is 
itself a leading source of disturbing force acting on the 
TM. Achieving the required extremely low level of spuri- 
ous acceleration requires identification and suppression of 
all interactions that can compete with gravity in defining 
the trajectory of the particle. The disturbances identi- 
fied thus far have been collected in an overall noise model 
[EHQthat has been used to optimize the LISA design. 

Particular attention has been dedicated to the interac- 
tion between the TM and the position sensor that sur- 
rounds the TM and is used to guide precision thrusters 
that keep the spacecraft centered around the free-falling 
TMs. As such, the design of the TM and position sensor 
- referred to together as the "gravitational reference sen- 
sor" or GRS - is critical to achieving the LISA force noise 
goals. The current GRS is a 46 mm cubic TM of Au-Pt 
(2 kg) surrounded by an array of conducting electrode 
surfaces used for a capacitive position readout and force 
actuation scheme (see Fig [1] and The TM material 

is chosen for low magnetic susceptibility and residual mo- 
ment, to minimize interaction with magnetic field noise 
9]. To reduce short range stray electrostatic effects, the 
distance ( "gaps" ) between the TM and the electrodes are 
relatively large, 4 mm on the interferometry axis, and 
all conducting surfaces are Au-coated to provide elec- 
trostatic homogeneity. To limit temperature differences 
across the TM, and thus also the temperature-gradient 
related forces which are the subject of this article, the 
sensor housing is made of a high thermal conductivity 
composite structure of molybdenum and sapphire. 

The unprecedentedly small level of force noise needed 
for LISA requires measurement of the small disturbing 
forces in order to verify the feasibility of the gravita- 
tional wave sensitivity goals. Given the possibility of 
force noise sources that are not accurately modelled or 
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X- face sensing electrodes 



FIG. 1: Schematic of the six degree of freedom capacitive po- 
sition sensor for LISA Pathfinder and LISA. The sensor elec- 
trodes are obtained by Au-coating of an insulating substrate, 
mounted on an external metallic frame. At left, sensing elec- 
trodes are in orange, sensor bias injection electrode are red, 
and the supporting housing is brown. The black circles repre- 
sent holes to allow access for laser hght and for a TM caging 
mechanism during launch. The x axis represents the sensitive 
LISA interferometry axis relevant for force noise. A internal 
section of the GRS is shown, with the blue arrows represent- 
ing the path linking the two X faces of the TM for gas and 
thermal photons. The TM-electrode gaps on the three axes 
are = 4 mm, dy = 2.9 mm and dz = 3.5 mm. 



perhaps not even previously identified, these tests should, 
as much as possible, use the final flight hardware config- 
uration, reflecting not only the nominal design materials 
and geometry, but also the relevant machining, surface 
finishing, impurities, and cleanliness. Such testing is be- 
ing pursued both in space, where the LISA Pathfinder 
mission will perform a dedicated flight test of the full 
LISA drag- free control system [l^, [ll|, and on ground, 
where torsion pendulum dynamometers can characterize 
the force noise generated inside the GRS [l^. A recent 
torsion pendulum study of the GRS capacitive sensor to 
be flown aboard LISA Pathfinder has demonstrated the 
absence of unknown surface forces to within roughly an 
order of magnitude of the LISA goal at 3 mHz [l^ . Other 
studies have measured important electrostatic and elastic 
coupling effects [13, [l^ [la [13, [l3| ■ Ongoing work aims at 
extending these studies to lower force noise levels, lower 
frequencies, and to even more representative hardware 
configurations. 

This paper describes a detailed study of the specific 
class of force noise induced by thermal gradients. Ther- 
mal gradient-induced forces have been modeled as part of 
the acceleration noise budget for LISA [i[i,[3,[il[23 and 
have been the subject of preliminary experimental inves- 
tigations [16, 17, ,21]. A temperature difference between 
the sensor surfaces on opposing sides of the TM con- 
verts into a net force through at least three mechanisms: 
the radiometric effect, differential radiation pressure, and 
temperature dependent outgassing. Broadly speaking, 
the radiometric and radiation pressure forces arise in the 



temperature dependence of the momentum transferred 
to a TM face by impacts of, respectively, residual gas 
molecules and thermal photons. If, at the LISA operating 
temperature (293 K) and pressure (10~^ Pa), we apply a 
simplified "infinite plate" model - to be discussed, along 
with its limitations and a more accurate simulation, in 
Sec. In] - the radiometric and radiation pressure effects 
can be estimated to produce a net force per degree 
of temperature difference, of roughly 20 and 30 pN/K, re- 
spectively, always pointing from the hot side of the sensor 
towards the cold. A third effect, more difficult to model 
and dependent on surface cleanliness, arises in the tem- 
perature dependence of outgassing from the sensor sur- 
faces and the net momentum thus imparted on the TM 
by desorbed molecules in the presence of a thermal gra- 
dient. As will be discussed in the following section, this 
effect has been expected to contribute roughly as much 
as the other two effects. The three effects respond coher- 
ently to the same temperature gradient and sum to give 
a net temperature-difference response of roughly « 
100 pN/K. 

This estimated "transfer function" sets the tolerable 
level of temperature difference fluctuations. In order to 
hold the thermal gradient-induced force noise to roughly 
10% of the overall force noise budget, the fluctuations 
in temperature difference across the GRS must be less 
than 10~^ K/Hz^/^. This level of temperature difference 
stability, which motivated the use of high thermal con- 
ductivity materials in the sensor, looks feasible but is 
challenging at the lowest frequencies, where the passive 
thermal filtering of solar radiation intensity fluctuations 
becomes progressively less effect ive[22l. [23j . 

The article aims at a characterization, through exper- 
iment and simulation, of the thermal gradient transfer 
function relevant to LISA free-fall. This study is 
important to free-fall for LISA and other precision force 
measurements for several reasons. First, while the rough 
estimates discussed above for the radiation pressure and 
radiometric effects are readily obtained from simple for- 
mulas, they are based on an infinite plate model, which 
is not particularly accurate for the geometry of the pro- 
posed GRS, where the dimensions of the TM are only 
roughly 10 times larger than the TM-sensor gaps. As 
will be seen from the results of the simulations presented 
in App. [X] and App. [B] these corrections, based on the 
geometry and, in the case of radiation pressure, surface 
reflectivity properties of the sensor, have a quantitative 
impact on the estimates of these noise sources for LISA 
and on the measurements presented here. 

An even more important motivation for an experimen- 
tal investigation lies in the large uncertainties with which 
we can estimate the outgassing effect, which will depend 
not only on the sensor geometry, but also on the amount, 
geometrical distribution, chemical nature, and history of 
the surface contamination inside the sensor. The amount 
to which the purity of free-fall is threatened by outgassing 
can only be reliably estimated through measurement with 
representative sensor hardware. The same argument can 
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be applied to any other possible unmodelled or unidenti- 
fied thermal related disturbance. 



It should be noted that knowing the impact of ther- 
mal gradients can have an impact on the overall LISA 
experiment design, particularly on the amount of ther- 
mal filtering needed, or even whether active temperature 
control is necessary. Additionally, the measurement of 
can be repeated in-flight during the LISA and LISA 
Pathfinder missions, a "calibration" that, combined with 
the appropriate thermometer readings, will allow a cor- 
relation analysis of the noise and even a subtraction of 
the disturbance from the gravitational wave time series. 



The experimental investigation was performed by 
means of a torsion pendulum, with a thin fiber suspend- 
ing a hollow LISA-like test mass on-axis with its center 
of mass, inside a GRS prototype. Two prototype sen- 
sors were tested, the current LISA and LISA Pathfinder 
design showed in Figll] with 4, 3.5, and 2.9 mm gaps, 
and an earlier prototype with smaller gaps, 2 mm on all 
axes (these are described in Sec. IIII A[) . The geometry 
of the pendulum is sensitive only to torques rather than 
the net forces that are relevant to along the LISA 
interferometer x axis. As such, the basic measurement 
technique, described in detail in Sec. Illli will be to ap- 
ply a "rotational temperature gradient" such as to excite 
a measurable torque that is proportional to the applied 
heat and, through modelling of the different effects, to 
the transfer function -j^- The measurements, of excited 
torque as a function of the temperature gradient pattern 
and presented in Sec. HVl are performed for various sen- 
sor average temperatures and for a range of pressures, 
which allows a separation of the radiometric and radia- 
tion pressure effects from their known functional depen- 
dencies. 



The pendulum configuration and the measurement of 
a torque (rather than a force) signal necessarily compli- 
cates the experimental analysis, demanding an interpola- 
tion of the relevant sensor temperature distribution from 
a limited array of thermometers and comparison with 
the torques calculated from the radiometric and radia- 
tion pressure models. In addition to complicating the 
analysis, the torque measurement imposed by the pen- 
dulum configuration means that the measurement will 
be blind to potential outgassing effects that act directly 
on the center of the TM faces, as there will be no "effec- 
tive armlength" to convert such forces into measurable 
torques. However, in spite of these analysis difficulties, 
discussed in Sec. IIII[ and the limitations of the torque 
measurement, to be discussed along with the study re- 
sults in Sec. |Vl relevant conclusions can be drawn con- 
cerning the role of all relevant contributions to the overall 
thermal gradient effect for LISA free-fall. 



II. MODELING OF THERMAL 
GRADIENT-INDUCED FORCES IN THE LISA 
GRS 

In this section we describe the three main physical 
processes that convert fluctuating temperature gradients 
into noisy forces on the LISA TM. For the radiometric 
and radiation pressure effects, we present both the simpli- 
fled models mentioned in the introduction and the results 
from numerical simulations using realistic sensor dimen- 
sions and surface properties. A thorough explanation of 
the simulations is given in Appendixes [X] and |B] 

In the discussion that follows, we assume the average 
GRS temperature to be Tq, with T(r) the temperature 
distribution on the internal surface of the position sen- 
sor, facing the TM at a position r. Given the high ther- 
mal conductivity of the sensor materials and the small 
expected thermal loads, we consider small temperature 
differences (5T(r) = (r(r) - Tq) < Tq. Additionally, as 
the thermal conductance across the TM is much greater 
than that between TM and sensor, at least at a pres- 
sure of 10~^ Pa, the calculations assume the TM to be 
isothermal at Tq. Finally, to connect the results with the 
transfer function relevant to LISA and discussed 

in the introduction, we define the average temperature 
difference ISTx = — T^^ ^ between the average tem- 
perature of the two opposing inner sensor surfaces or- 
thogonal to the X axis, with the positive and negative X 
surfaces defined X+ and X_ as in Fig[TJ While for LISA 
only the force component along the sensitive x axis, F^, 
is relevant, here we also calculate the z component of 
the torque, Nz, that coincides with the torsion pendu- 
lum fiber axis and is thus relevant to the measurements 
presented in this article. 



A. Radiometer effect 

In gas systems where the mean free path is long com- 
pared to the container dimensions - such as for the LISA 
GRS at 10~^ Pa - the hydrostatic equilibrium condition 
of pressure uniformity, which is valid in the dense gas 
limit, is no longer relevant. In the presence of a temper- 
ature gradient, with the equilibrium condition that there 
be no net flux of molecules between hot and cold zones, 
a stable pressure difference is established. The radiome- 
ter effect refers to this steady state pressure difference 
that, for LISA, arises in the net difference in momentum 
transferred to opposing faces of the TM in the presence 
of a temperature gradient. 

Radiometric forces can be estimated with the tran- 
spiration theory for a gas in the free-molecular regime 
0) IHj HI] J which considers that the momentum distri- 



^ In this sign convention, ATx positive will create a positive Fx, 
and thus the various contributions to ^."^^ ^ will be positive. 
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bution of molecules arriving on a surface element depends 
on the temperature of the last surface from which they 
departed, and will thus be affected both by the system's 
specific temperature distribution and geometry. We con- 
sider here the pressure Ppi between two parallel plates, 
"infinite" in the sense that their separation is negligible 
compared with their linear extent. With the two plates 
at different uniform temperature Ti and T2 and in equi- 
librium with a particle reservoir at pressure P and tem- 
perature Tq, the approximate solution is 



(1) 



This result holds even in the case that the molecules are 
not fully "thermalized" upon colliding with the surfaces 
and is only slightly modified if the effective thcrmaliza- 
tion factor is different between the two surfaces 25]. In 
the limit of TM-sensor gaps that are much smaller than 
the effective spatial extent of the temperature perturba- 
tions, we apply this formula to obtain the force dFfj{r) 
( "R" for radiometric) which acts normally on TM surface 
element ds at position r : 





Pds- 



P 6T{r) 
T To 



ds (2) 



Here we have used the above mentioned approximation 
of an isothermal TM and ST{r) < Tq. This simphfied 
"infinite plate" model assumes that all (and only) the 
molecules emitted by the GRS surfaces directly in front 
of the TM (for the X faces, these are the red and yellow 
zones in Fig[2]) strike the opposing TM face. As a conse- 
quence, edge effects are negligible, and the net shear force 
acting parallel to any TM surface is zero. Thus, the in- 
tegrated radiometric contribution to the x force comes 
from the GRS surfaces on the and X_ faces and is 
given by 




{T)5T{y,z)dydz 



(3) 



where, as defined above, Tx+ — = is the aver- 
age temperature difference between the X+ and X- GRS 
surfaces, and A is the the area of a TM face. For exam- 
ple, in the idealized temperature difference distribution 
shown at left in FiglSl AT^ = 25T. 

The radiometric "transfer function" is thus 



dF,, 



d{AT,) 



AP 
18pN/K X KR 



P \ /293K 



10-5 Pa/ V 7c 



(4) 



Here, we include a radiometric correction factor kr in 
order to incorporate modifications to this simple model 
that arise when considering a realistic sensor geometry 
with finite dimensions, which will be discussed shortly. 
The factor kr will be a function of the sensor geometry, 
but also of the particular temperature distribution inside 
the sensor that produces a given average temperature 
difference AT^. 

In this model, the torque component Nz, relevant to 
our torsion pendulum experiment, is calculatd by inte- 
grating the moment of the normal radiometric force el- 
ement in Eqnl2] over the area of the x and y TM faces, 
yielding 



(±) ST{y, z)ydydz 




/ (t) 5T[x, z)xdydz 
■ V JA 



(5) 



The sum of the integrals within the parentheses 
will be called "thermal integral" and indicated with 
Jg T{s)b{s)ds. The arm-length b{s) coincides, taking the 
appropriate sign, with the coordinate y on the X faces 
and the coordinate x and on Y faces. Thus Eqn[5] can 
be simply written as an integral over the TM X and Y 
surfaces, 



Nz c 7fl-^ / T{s)b{s)ds. 

4Jo JS 



(6) 



By analogy with the correction factor kr introduced in 
Eqn|4l we include here the factor 7^^, as a radiometric 
torque correction factor for a finite size sensor. 7^ will 
also be a function of the sensor geometry and of the par- 
ticular temperature distribution inside the sensor that 
produces a given thermal integral. 

To illustrate the significance of the "thermal integral," 
we consider the simple temperature distribution on the 
right in Figl2l which is an idealization of the temperature 
pattern imposed in our experiment. In this case. 



T{s)b{s)ds = ST^^^iST^^. (7) 



Each "electrode" (the half surfaces of the X faces) con- 
tributes with a torque proportional to the temperature 
difference ST multiplied by its area L x -j and the " arm" 
■J (that is the distance of the center of the electrode from 
the center of the TM face). For the same ST, the thermal 
integral scales as the cube of the linear size of the TM. 

For realistic sensor dimensions, finite size effects begin 
to erode the infinite plate assumptions used here. First, 
not all molecules leaving the GRS surface directly in front 
of the TM (for the X faces, the yellow and red zones in 
Figl2) strike the TM. Other molecules originating in the 
corners (the grey corners in FiglJ) do strike the TM, on 
the X faces but also on the Y and Z faces. Additionally, 
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FIG. 2: Idealized temperature gradient configurations. The 
orange square represents a section of a cubic TM of size L at 
uniform temperature Tq. Yellow and red parts represent do- 
mains of uniform temperature on the inner surface of the elec- 
trode housing, as indicated next to them. Left panel: because 
of radiometer effects (and, as shown in Sec. Ill Bl radiation 
pressure), the illustrated temperature pattern is expected to 
induce a net force Fx- Right panel: the illustrated tempera- 
ture pattern is expected to induce a torque around the TM 
z axis. In the simplified infinite plate model of the radiometer 
and radiation pressure effects, the gray parts of the electrode 
housing do not contribute to F^, even in the presence of a 
temperature gradient (see Eqn[3}. The gray part can, in this 
simplified model, contribute to A''^ , in the event of a thermal 
gradient, through the forces that they cause orthogonally to 
the TM faces in front of them (as in EqnO. 



molecules striking the Y and Z faces impart shear forces 
along X, which in the presence of a temperature gradi- 
ent is rendered asymmetric by the finite TM dimension 
(illustrated in Fig 13]). These phenomena are enhanced 
with larger gaps and, in connection with temperature- 
induced asymmetry of the molecular velocity distribu- 
tion, are expected to contribute to both and iVfl^. 
To understand the role of these corrections to the simple 
radiometric model presented above, we have performed 
a numerical gas dynamics simulation, which is described 
in detail in App. [X] 

Simulations have been performed for a cubic 46 mm 
TM inside a rectangular box, with several different gap 
sizes, including the geometry for the LISA Pathfinder 
sensor design (shown in Fig[T]). The simulations con- 
sider a temperature difference IST^ constituted by 
and X- faces that are at constant temperature (differ- 
ing by ATj), with a linear temperature gradient along 
the Y and Z faces (as illustrated in Figla]). For the TM 
force along x, the simulations indeed verify the simplified 
model of Eqn|3](Kfl = I), to within 10%, for the smallest 
TM-GRS gaps studied (1 mm, see Fig [H] in App. E|, 




FIG. 3: Cartoon of the GRS cross section, illustrating shear 
forces along x caused by molecules striking the Y TM faces 
in the presence of a linear thermal gradient in x (as in the 
simulations of App. Molecules striking the Y faces im- 
part shear forces parallel to the TM surface along x. For a 
non-infinite TM, some molecules (such as those represented 
by squares) emitted from a given point on Y sensor surfaces 
do not strike the TM. As such, this asymmetrizes the x mo- 
mentum imparted to the TM by the molecules originating 
in the given sensor-face location, which, in the presence of 
a thermal gradient produces a net thermal shear force (blue 
arrows). This is reinforced by molecules (represented by tri- 
angles) starting from the X faces near the corners. From the 
standpoint of a surface element of the TM, there is a net im- 
balance in the shear momentum imparted by the molecules 
arriving from the hot (red) and cold (yellow) sides that con- 
tributes to the overall translational radiometric effect. This 
shear effect increases with increasing gap size. 
Note that the momentum components transferred normal and 
parallel to the Y surface contribute to the torque around the 
z axis with opposite sign. As such, these shear forces are also 
responsible for a reduction in the radiometric torque, relevant 
to the experiments presented here. 



with the force coming nearly exclusively from molecules 
impacting the X faces. With increasing gap size, how- 
ever, the calculation demonstrates the role of shear forces 
along the Y and Z faces, which contribute to an excess 
of roughly 40 % above the infinite plate model predic- 
tion for the largest gaps of the study (7 mm). For the 
LISA Pathfinder GRS design, the radiometric force trans- 



fer function 



d(AT) 



exceeds the infinite plate prediction 



EqnlUby roughly 25 % ( /v_r « 1.25), with a statistical 
error of several percent (results are summarized in Ta- 
ble [H along with those from a similar radiation pressure 
simulation) . 

The simulations employing a linear temperature gra- 
dient along X also allow an estimate of the radiometric 
torque relevant to the torsion pendulum experiment, as 
the linear temperature profile along y produces a torque 
around the z axis. Here, a net decrease of the radiometric 
effect is observed with respect to the simplified model of 
Eqnl5](see Fig[Tn]of App. E|. With increasing gap, there 
is a suppression of the radiometric torque, dominated by 
the same shear forces mentioned above, which tend to 
work against the torque contribution of the dominant 
normal component of the molecular impacts (illustrated 
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1.25 




a = 1 


1.17 


K.RP 


a = 0.05 (specular) 


0.32 




a = 0.05 (diflFuse) 


0.75 



TABLE I; Results of the thermal gradient simulations, where 
kr and krd represent correction factors to the infinite plate 
model for, respectively, the radiometric effect (Eqn|4ll and dif- 
ferential radiation pressure (Eqn llO|l . The simulations assume 
the LTP GRS design geometry and a temperature difference 
between two opposing GRS X faces, each at a uniform tem- 
perature and with a linear temperature gradient along the 
Y and Z faces. Results for radiation pressure are shown for 
different values of the absorption a and nature of reflection 
(specular or diffuse). Statistical errors are of order several 
percent. 



in Figl3]). The magnitude of the torque reduction, with 
respect to the infinite plate model, is in fact roughly con- 
sistent with the shear force -relevant to the increase in 
radiometric force reported above -multiplied by an arm 
length of L/2, or half the TM width. For the Y face 
torque contribution for the 2 mm gap sensor, we find 
roughly 83% of the torque calculated with the infinite 
plate model (EqnlHl), and only 65% for the larger-gapped 
LTP GRS prototype. 

The true correction factor will depend on the real 
temperature profile for the experiment, rather than the 
contribution of a single face that sees a linear temper- 
ature profile. However, the experimental uncertainty in 
estimating the true temperature profile in these experi- 
ments (see Sec. IIII CI and App. [C]), particularly in the 
corners of the sensor where the finite-size corrections are 
most relevant, limits the value of a simulation performed 
with a "full" (but very approximate) sensor temperature 
distribution. A linear temperature distribution also turns 
out to be a decent approximation for the experimental 
temperature distributions. As such, the torque results 
obtained for a linear temperature gradient, extrapolated 
to a full sensor, yield indicative estimates, ^b. ~ 0.65 and 
7fl ~ 0.8, for comparing the experimental results for the 
torque measurements for, respectively, the LTP GRS and 
2 mm sensor prototypes. 



B. Thermal Radiation Pressure 

The effect of thermal radiation pressure can be evalu- 
ated from the momentum transfer of thermal photons 
emitted by all radiating surfaces inside the GRS. We 
again assume an infinite parallel plate approximation, 
in which all thermal photons wind up being absorbed 
in the vicinity of emission, either on the parallel surface 
opposite the point of emission or reflected back onto the 
original emitting surface. Assuming the same emissivity 
for all the involved surfaces, each surface element ds on 
the internal TM-facing surfaces of the sensor contributes 



with a force dFupij) acting on the TM in the direction 
perpendicular to the surface element itself. This force is 
given by 

dFnp{v) ^ YT\v)ds ^ §r4(r)ds + ^T3<5r(r)ds(8) 

with cr the Stefan-Boltzmann constant and c the speed 
of light. Assuming that the TM is isothermal, and thus 
contributes no net recoil force from its radiation, the net 
force acting on the TM in the sensitive x direction for 
LISA will be proportional to the average temperature 
difference between opposing sensor faces, given by 



(9) 



where Eqn[5] is integrated over the EH surfaces that di- 
rectly face the TM (the yellow and red zones in the left 
image m Figl2|)i. Thus, the radiation pressure transfer 
function can be expressed 



dF, 



diAT,) 



RP 



..pa|t3 



27pN/K X Kp^p 



To 
293 K 



(10) 



As in Sec. IH Al we can integrate the moment of the 
force element dFpp (Eqn[5]) to calculate the torque N^, 



Nrp.~-Irp^T'^ I T{s)b{s)ds 



(11) 



where both the X and Y faces contribute. We note that 
in analogy with the radiometric effect in the previous 
section, we have introduced, in Eqns. 1101 and [TTl respec- 
tively, we have introduced the radiation pressure force 
and torque factors Kpp and 7_rp, to account for correc- 
tions due to the inaccuracy in applying this simple infi- 
nite plate model to our true sensor geometry. 

As with the radiometric effect, the assumptions for this 
simple radiation pressure model begin to break down 
when the gap sizes are not negligible compared to the 
TM dimensions. Just as for gas molecules, there are 
corrections due to emitted photons "missing" the TM 
and hitting adjacent faces, as well as adsorbed photons 
along the lateral {Y and Z) faces imparting a shear 
force and associated torque on the TM. In addition, the 
presence of finite reflectivity (absorptivity a < 1 means 
that bouncing thermal photons will hit multiple surfaces, 
imparting force and torque components with different 
signs, before adsorbing. This latter effect is particu- 
larly significant for the metallic surfaces inside the sensor. 



^ Eqn|9]is independent of the absorptivity coefficient if the TM and 
sensor wall absorptivities, aj^M and ag respectively, are equal. 
In the event that they are not equal, Eqn[2]is modified by a factor 



°s(2 — utm) 
'^TM+'^S—"-TM"'S 



which can range from to 2. 
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which are expected to have a high 90-95 % reflectivities 
(a ~ 0.05 — 0.1) at thermal photon wavelengths [l^], and 
systematically serves to homogenize the radiation pres- 
sure felt by the TM, suppressing the differential radiation 
pressure effect. 

To properly analyze such effects, a numerical simu- 
lation was performed and is described in detail in Sec. 
IB] The simulation studied several representative sim- 
plified sensor geometries and sensor temperature profiles 
representative for the force for LISA and for the torque 
relevant to the experiments presented in the following 
sections. Additionally, different extreme cases for the in- 
frared absorption a and for the nature of the refiection 
(specular or diffuse) were considered. For the force, a 
simple linear temperature gradient along the x axis was 
assumed, with uniform temperature across each X face. 
In the limit of vanishing TM-sensor gaps (0.2 mm), the 
infinite parallel plate model for the force (EqnlQ]) was con- 
firmed, with the factor |^ (or kjip « 1), to within several 
percent, even for absorption as small as 0.1. For perfect 
absorption (a = 1), there is an increase in the radiation 
pressure effect akin to that for the radiometric effect, 
with Knp = 1.17 for the LISA Pathfinder geometry (see 
the summary of results in Table |T| . For the case of the 
high reflectivity surfaces expected for the LISA sensor - 
using a = 0.05, with specular reflection - there is a sup- 
pression of the effect, with k^p — 0.32. This result thus 
indicates a significant, if not ordcr-of-magnitude, change 
to the estimate of the effect given in simplified radiation 
pressure models for LISA [5, 6, J,]. 

For the torque, assuming temperature profiles similar 
to that estimated in our experiments (using the technique 
described in Appendix A), the suppression of the differ- 
ential radiation pressure effect is more pronounced, with 
photons moving, in only a couple bounces, from one half 
of the TM to the other (or to an adjacent face), inverting 
the sign of the resulting torque. For the high refiectiv- 
ity surfaces (a = 0.05 and specular refiection) reasonable 
for the gold coated TM and sensors used, we obtain ra- 
diation pressure torque correction factors jpp of order 
0.4 and 0.3 for the 2 mm and 4 mm (LISA Pathfinder 
design) prototypes used in the experiment. This signifi- 
cant reduction has an impact on the interpretation of the 
experimental results discussed in Sec. IIVI 



C. Asymmetric temperature dependent outgassing 

Outgassing of molecules absorbed by the internal walls 
of the GRS increases the residual pressure surrounding 
the TM. Additionally, an asymmetry in the molecular 
outflow can be created or enhanced by the temperature 
dependence of the outgassing, resulting in a differential 
pressure that is driven by the GRS temperature gradient 
itself. 

As asymmetric outgassing depends on the amount and 
type of absorbed impurities inside the GRS, it is a con- 
tamination effect that does not lend itself easily to a 



flrst principles calculation, in contrast with the previ- 
ous two effects. However, the gas flow from a sur- 
face can be modelled with a temperature activation law 
Q{T) — Qoe*-~®^"^\ where Qo is a flow prefactor and 
O the activation ternperature of the molecular species 
under consideration [28l [29| . Our simplified model con- 
siders only a single outgassing species which is uni- 
formly desorbed from identical GRS surfaces, and that 
the molecules emitted from the sensor walls collide with 
the TM without sticking. A temperature gradient AT^ 
would induce an asymmetric rate of outgassing across the 
TM AQ(T) w Q{To) (O/To) (ATjTo). The consequent 
difference in pressure APout across the TM would cause 
a force Fqgx given by 



F, 



OGx 



AAPoG 



^AOm^^Q^^^ (12) 



C, 



eff 



where Cef / is a geometrical factor resulting from a com- 
bination of the conductance of the paths around the TM 
and through the holes in the GRS electrode housing walls, 
estimated to be roughly 4.3 xl0~^ m^/s for the LISA 
/ LISA Pathfinder GRS geometry Q. The asymmetric 
outgassing transfer function is then given by 



dAT, 



A 



OG 



40 



QjTo) e 

pN f QjTo) 
K Vl-4nJ/s 



e \ /293KV 
3xl04Kj 



We note that the temperature dependence of the effect, 
partially hidden in (5(To), is ^exp(— 0/T), at least in 
this simplified model with a single outgassing species, 
and is thus a rapidly increasing function at the temper- 
atures of our study. The values used for the flow Q(To) 
and the activation temperature correspond to typical 
room temperature behaviors for the metal and ceramic 
materials used in the GRS [s^l . This source of cou- 
pling to the thermal gradient thus appears to be roughly 
equal in magnitude to that of the combined radiometric 
and radiation pressure effects (Eqns. Hland fTU)) . However, 
given the likely large variability due to materials, assem- 
bly, cleaning, and thermal history of the sensor, this noise 
source has a much larger uncertainty and is a key moti- 
vation for measurement of the overall thermal gradient 
transfer function. If the outgassing phenomenon were in- 
deed uniform across the GRS surfaces, the torque transfer 
coefficient could be found by integrating the moment due 
to the differential outgassing pressure over the sensor sur- 
faces, as for the two preceding thermal gradient effects. 
This would allow conversion between force and torque 
through an effective arm length, of order ^ for a simpli- 
fied temperature distribution like that in Fig[2l However, 
unlike the previous two effects, the numerous features in 
the sensor - holes, interfaces between the metal and ce- 
ramic surfaces - suggest that the outgassing could be 
concentrated in specific localized areas, rather than uni- 
form. As such, there is no simple, unequivocal extrapola- 



tion of 



dAT^ 



OG 



from a torque measurement. However, 



the torque measurements are sensitive to nearly all pos- 
sible outgassing effects, regardless of their model, and 
the lack of a significant excess beyond the well-modelled 
radiometric and radiation pressure effects will certainly 
build confidence that there is no dominant asymmetric 
outgassing effect for LISA. 



III. EXPERIMENTAL DETAILS 
A. Gravitational Reference Sensor prototypes 

This experimental campaign studied two GRS proto- 
types realized as part of the sensor design development 
for LISA. They feature different electrodes geometries, 
construction techniques, and TM sizes. Both sensors are 
based on a Gold coated Molybdenum housing, with in- 
sulating elements made from the Aluminum-based, high 
thermal conductivity ceramic SHAPAL. The first proto- 
type sensor, referred to here as "TN2mm" uses a 40mm 
cubic TM, with a gap - the distance between the TM 
and the surrounding electrode and electrode housing sur- 
faces - of 2 mm, symmetric on all three axes. The elec- 
trodes are Au-coated bulk Mo plates separated from the 
grounded housing by bulk Shapal plates. The second 
prototype, "EM4mm," reflects the final GRS design for 
the LTP mission, with a 46mm cubic TM and gaps of 4, 
2.9 and 3.5 mm on, respectively, the x, y, and z axes 
(see FigH]). In contrast with TN2mm, the electrodes for 
EM4mm are obtained by Au sputtering deposition di- 
rectly on the bulk insulating SHAPAL plates. 
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FIG. 4: Top panel: prototype TN2mm was equipped with five 
thermometers and the experiment was performed by alterna- 
tively switching on and off the heaters HIE and HIW (red 
circles). Bottom panel: prototype EM4mm, equipped with 
eight thermometers. For prototype EM4mm we employed 
all four heaters, alternatively switching the pairs (H1E,H2W) 
and (H2E,H1W). We note that T/afce is included here on the 
TN2mm X-\- face to indicate the position of a temperature 
reading that, as described in App. [C] has been artificially 
added to help estimate the temperature profile, in the ab- 
sence of a thermometer near to heater H2W. 



1. Electrical heaters and thermometers 

Electrical heaters have been installed on the exter- 
nal surfaces of both sensor prototypes, in correspondence 
with the X face electrodes. They are displaced with re- 
spect to the pendulum torsion axis (see Fig|3]) and can 
thus apply the rotational temperature differences that 
convert the forces described in Sec. [ll] into torques, to 
which this pendulum is sensitive. The heaters are made 
from of twisted pairs of Manganin wires several meters 
long, wound onto 5 mm diameter Al cylinders. The 
heaters are attached to the sensor with a high thermal 
conductivity, vacuum compatible glue, for the TN2mm 
prototype, and tightly screwed onto the GRS surfaces 
for EM4mm. The resistance of each heater is 95 ± 1 
Ohms. A 1 Hz square wave voltage, up to a maximum 
of 7 V across the resistor (nearly 0.5 W), is supplied to 
the heaters with high output current operational ampli- 
fiers. Twisting the heater wires limits the magnetic field 
produced by the heater current, and the use of the 1 Hz 
square wave allows application of constant heater power, 
while at the same time converting any residual magnetic 
field-related torque to a frequency well above our mea- 
surement frequency of 0.4-0.5 mHz (see Sec. IIII CI) . 



With the aim of reconstructing the temperature pro- 
file of the GRS, temperatures at different positions are 
measured for each sensor, by a set of PtlOO thermome- 
ters on the electrode housing external surfaces. They are 
glued to the sensor by means of a high thermal conduc- 
tivity vacuum compatible glue. Five thermometers were 
used on the TN2mm prototype and eight on EM4mm. A 
sketch of the thermometer and heater locations is shown 
for both sensors in Fig[4l 



B. Torsion pendulum apparatus 

This experimental investigation has been performed 
with a torsion pendulum facility [12, 0, Q El [3 IH , 
where a representative copy of the LISA test mass is sus- 
pended by a thin fiber and hangs inside a prototype of 
the GRS. The small torsional elastic constant of the fiber 
allows for high resolution measurement of torques around 
the fiber axis, via measurement of the pendulum angular 
deflection. Measurements of the pendulum angular de- 
flection noise can be used to study the force noise acting 
on the LISA TM. Additionally, controUedly modulation 
of known disturbance fields, such as the thermal gradi- 
ents in this study, with coherent detection of the resulting 
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FIG. 5: Schematic of the torsion pendulum apparatus. The 
TM is the main part of the inertial member of a torsion pen- 
dulum hanging from a long torsion fiber designed to have the 
minimum stiffness and thus the maximum sensitivity. The 
electrode housing surrounds the TM and is attached to a pre- 
cision positioner that can center it relative to the TM in 6 de- 
grees of freedom. 



pendulum deflection, can allow a study of known force 
noise sources relevant to LISA free-fall. 

The pendulum is designed for minimum elastic stiffness 
and thus maximum torque sensitivity. A light-weight hol- 
low TM (Au-coated Ti for TN2mm, and Au-coated Al for 
EM4mm) has been chosen to maximize the torque sensi- 
tivity, allowing use of a thin (25 /xm) tungsten fiber, with 
a resulting resonant frequency near 2 mHz and quality 
factor near 3000. The pendulum and sensor are mounted 
on independent micro-positioners, allowing a 6 degree of 
freedom adjustment of the relative position of TM inside 
the electrode housing. The torsion pendulum deflection 
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FIG. 6: Illustration of the thermal gradient effects measure- 
ment technique for the the TN2mm prototype. Here 200 mW 
are alternatively applied to heaters (HIE and HIW in Fig|4]) 
located off-axis on opposing sensor faces, switching at 1000 s 
intervals (top panel). The central panel shows the induced 
temperature difference AT = (T2 — T4). The resulting ther- 
mally induced torque causes the pendulum to respond coher- 
ently, as shown in bottom panel, with an angular deflection 
at the heater modulation frequency superimposed on the pen- 
dulum free oscillation. 



angle was read out, with 10 Hz sampling rate, by an op- 
tical autocollimator measurement of a mirror that was 
part of the torsion member. 

The vacuum chamber that accomodates the pendulum 
and the GRS is evacuated by means of a turbomolecu- 
lar pump, and its pressure is monitored by means of an 
ionization gauge (Varian UHV 24) . The whole apparatus 
sits in a thermally controlled room with 50 mK long term 
stability. 

The pendulum torque resolution relevant to these mea- 
surements is typically 5-10 fN m / Hz^/^ near the 0.4 - 
0.5 mHz heater power modulation frequency, which is a 
factor 2-3 above the Brownian thermal noise limit for the 
pendulum. For the several hour integration times used in 
these measurements, this equates to a torque resolution 
of order 0.1 fN m. 



C. Measurement technique and data analysis 

The general experimental technique, illustrated in 
FiglH is to create an oscillating asymmetric tempera- 
ture distribution by modulated heater powers and then 
measure the resulting coherent oscillations both in the 
torque on the TM and in the GRS temperature dis- 
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tribution. Varying and monitoring the average sensor 
temperature and pressure, in addition to the oscillating 
temperature difference pattern, allows separation of the 
different thermally induced torque effects for their de- 
pendence on these parameters. In presenting the mea- 
sured torque data, they will be normalized for the ampli- 
tude of the applied temperature gradients as expressed 
by the estimated applied "thermal integral," the integral 
JgT{s)b{s)ds introduced in Sec. to characterize the 
radiometric and radiation pressure torques in the infi- 
nite plate limit. This normalization allows comparison 
between measurements of different thermal gradient am- 
plitudes and comparison with the models for the radio- 
metric and radiation pressure effects discussed in Sec. HH 



used. However, given the coarse sampling of the temper- 
ature distribution with only 5-8 thermometers and the 
approximate nature of the analysis used, we have chosen 
not to "correct" the data presented for this systematic 
error, but rather to present the data obtained from the 
relatively straightforward interpolation. The uncertain- 
ties in the estimation of the thermal integral are of order 
20% and 50% for, respectively, the EM4mm and TN2mm 
sensors. Given the repeated use of the same heating pat- 
tern for each of the two sensor datasets, the uncertainty 
in the thermal integral estimation can be thought of a 
single scale factor for each dataset. 



2. Torque estimation 



1. Rotational thermal gradients: application and estimation 

For the TN2mm prototype, the oscillating temperature 
pattern were obtained by alternately switching on and off 
two opposing heaters for 1000 s, as illustrated in FigH] 
For the EM4mm prototype, two pairs of opposing heaters 
are switched on and off for 1250 s. These modulation fre- 
quencies, 0.5 and 0.4 mHz respectively, result in coher- 
ent pendulum torques that are below the roughly 2 mHz 
resonant frequency. Typical peak to peak temperature 
differences across the sensor ranged from tens to hun- 
dreds of mK. The alternate switching of the heaters pro- 
duce an oscillating temperature gradient pattern across 
the EH, with the time dependence of a lightly filtered 
square-wave. We note that the thermal time constant 
for equilibration across the GRS, relevant to the temper- 
ature difference (as plotted in Figd]) and to the thermal 
integral, is of order 100 s. 

The temperature pattern was recorded in the readings 
Ti{t) of each thermometer i at its position (si, 52)1, where 
si and S2 are the two coordinates mapping the 4 X and 
Y GRS surfaces. For each cycle k of the heater power 
modulation, we calculate the component T{si, S2,oj, k)i 
of Ti{t) at the modulation carrier frequency and at its 
third and fifth harmonics, where the filtered square wave 
thermal pattern is known to have significant spectral con- 
tent. An estimation of the components T(si, 52,^^, k) at 
a generic position (sl,s2) was then obtained by spatial 
interpolation of the values T{si,S2,uj, k)i (for details see 
the appendix App. [C]). For each demodulation cycle 
k, it was then possible to obtain an estimation of the to 
component of the "thermal integral" , as defined in Eqn[5] 
and 



T{s)b{s)ds 

S / (ui,k) 



The thermal gradient induced torque on the TM is 
detected in the coherent deflection of the pendulum. The 
torque components are calculated for each cycle fc, at the 
odd harmonics of the heater power modulation frequency, 
using the torsion pendulum transfer function, as in Refs. 

and [l3|. The resolution of order O.lfN m obtained 
for a several hour integration corresponds to a signal-to- 
noise ratio of order 10 for the smallest thermal gradient- 
induced torques measured. 

The torque data for each cycle are normalized by divi- 
sion by the calculated thermal integrals. 



N, 



(uj,k) 



(w,fe) 



{J,T{s)b{s)ds) 



(15) 



(ui.k) 



The units of N, 



{u,k) 



Pa 
K ■ 



3. Average temperature: control and estimation 

Average temperatures between 283 K and 314 K were 
obtained by changing the temperature of the room con- 
taining the pendulum apparatus and by changing the to- 
tal DC power applied to the heaters. The average tem- 
perature of the electrode housing and TM Tq is derived 
by averaging the thermometer data for each measure- 
ment cycle. In particular, average temperature data for 
UTN2mm were obtained with Tg = ^^-5^, whereas for 



EM4mm we used Tr 



4 



(see FigHl 



4-. Average pressure: control and estimation 



V J A 



The pressure in the vacuum enclosure were set to differ- 
ent equilibrium levels, from roughly 2 x 10~^ Pa to lO""* 
/ (ip) T{x z uj k)xdydz (iJ^^' changing the conductance between the vacuum 
J A ' ' ' vessel and the pumping system. The pressure in the vac- 

uum vessel was measured by an ionization gauge (Varian 
UHV 24) , sampled at 1 Hz and averaged for each cycle of 
the heater modulation. The calibration of the ion gauge, 
for which no certification was available, was obtained by 



As discussed in appendix App. [Cl our analysis shows that 
the simple interpolation technique used underestimates 
the thermal integral, by 20-30 % in the simplified analysis 
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testing against a similar device (PfeifFer UHV) , that was 
certified by the manufacturer to within 15%. It was not 
possible to use the Pfeiffer gauge constantly because of 
its very high filament current, which produced an intol- 
erable electrostatic charging of the TM. The calibration 
was performed by simultaneously acquiring both gauges 
while slowly changing the pressure in the vacuum vessel 
with a valve on the pumping line, over a pressure range 
from roughly 5 x 10~^ Pa up to 2 x 10"** Pa, roughly the 
same range spanned by our study. The resulting fit of 
the obtained curve was applied to calibrate the Varian 
gauge in all measurements. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

Results of the experimental investigation for both the 
TN2mm and EM4mm sensors are shown in Figl?! where 
we report N*^ ^-j , the component of the normalized torque 
at the heater modulation carrier frequency w, as a func- 
tion of pressure and for several average sensor temper- 
atures. The data are dominated by the linear pressure 
dependence of the radiometric effect, which will be dis- 
cussed in the following paragraphs. 

In FiglHwe plot the normalized torque data at the first, 
third, and fifth harmonics of the modulation frequency 
(0.4, 1.2 and 2.0 mHz) - each normalized for the thermal 
integral at the relevant frequency - with the 298.9 K data 
for the EM4mm sensor chosen as an example. The data 
at the three frequencies are equivalent within the sta- 
tistical errors of the measurement. In this study, no fre- 
quency dependence in the transfer function between mea- 
sured thermal gradients and resulting torque has been 
observed. Also, by varying the amplitude of the modu- 
lated heater power, we can check the linearity of the ther- 
mal induced torque effects in the applied thermal integral 
amplitudes. One example is shown in Fig|5]for the proto- 
type TN2mm, showing that the normalized torque data 
are reproduced, within the statistical errors, for the two 
different temperature gradient amplitudes. These two 
observations confirm the linear, frequency-independent 
nature of the thermal gradient torque models discussed 
in Sec. M 

In the following two sections we discuss the experimen- 
tal results in light of the models for the radiometric, ra- 
diation pressure, and temperature-dependent outgassing 
effects discussed in Sec. [Til From Eqns. [6] and [TTl we 
expect that 



N, 



P 



8ct, 



3c 



' 



(16) 



with an additional pressure independent term coming 
from temperature-dependent outgassing. By fitting the 
data for N*{Tq,P) in Fig[7] to the linear pressure de- 
pendence at each average GRS temperature Tq, we first 
study the radiometric effect, through the pressure slope. 
Finally through the zero-pressure intercept, we can study 
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FIG. 7: Normalized torque N*^ of Eqn IlSI as a function 
of pressure at different average temperatures. Top panel the 
data for TN2mm prototype sensor, the bottom panel the data 
for GRS Engineering Model for the LTP , indicated with 
EM4mm. Dashed lines are the results of linear least square 
fits to the pressure at each average temperature. Only data 
for the first harmonic of the excitation frequency are shown. 



the radiation pressure effect and any outgassing related 
phenomena. 



A. Comparison with the radiometer effect model 

According to EgnfTBl for the radiometer effect 
N*(u}, k) should depend linearly on pressure with a slope 
given by 



dP 



(17) 



where the correction factor 7/j = 1 in the limit of the in- 
finite plate model. Figl7]does indeed show that N*{lo, k) 
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FIG. 8: Normalized torques A'^ ^ as function of the residual 
pressure P, at the average temperature To « 298. 9K, calcu- 
lated at the fundamental frequency lu; — 0.4mHz (blue dots) 
and its odd harmonics 3a; — 1.2mHz (red) and 5a; = 2mH[z 
(black) (data for the EM4mm sensor). Linear fits also are 
shown for comparison. The three measured slopes and the 
three intercepts agree with each other within their statistical 
errors (smaller than 4%). 
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FIG. 9: Measurement of N^^ as a function of pressure, with 
prototype TN2mm, for different values of applied heat. Red 
points; alternating 200 ± 2 mW between heaters HIE and 
HIW. Blue points: apply 50 ± 1 mW continuously to both 
heaters and alternating 100 ± 1 mW. Lines are best fits to 
data. 



depends linearly on P, with a slope that decreases with 
To (barely visible by eye in the plotted data). 

In Fig[TOl we plot the slope ( ^gp' ) multiplied by 

V /To 

4:Tq, which thus represents the measured 7;?, a constant 
according to our model. The results yield 7^ = 1.34 
for the TN2mm prototype and 7^ — 0.92 for EM4mm. 
The RMS fit residuals are of order 1.5% for each sen- 
sor, resulting in overall uncertainties in 7/5 slightly below 
1%. We note that, for both sensors, the observed scatter 
in these data exceeds that estimated by the statistical 
torque uncertainty, by a factor 4-8. As such these data 
are not limited by the torque resolution, but likely by 
systematic uncertainties in the temperature integral and 
pressure measurements between points at different tem- 
perature. We note that the temperature ranges for both 
sets of data are sufficient to distinguish the temperature 
dependence, and the fit to the model (Eqn fTTl) is sig- 
nificantly better than a fit to a constant (temperature 
independent) or to 7^. 

Figure ]W\ also presents the values for jji calculated in 
the radiometric simulations for the two sensor geometries, 
assuming a simplified linear temperature gradient (Sec. 
El, which yielded = 0.83 for TN2mm and 7^ = 0.65 
for EM4mm. In comparing the experimentally measured 
values to these predictions, we have to consider the sys- 
tematic uncertainties in the estimate of the thermal in- 
tegral that has been used to normalize the torque data. 
The temperature interpolation scheme used to estimate 
the thermal integral, discussed in Appendix C, is likely to 
underestimate the thermal integral. By one approximate 



analysis, this underestimate could be roughly 20 — 30% 
for the two sensors. Underestimating the thermal integral 
would systematically increase the values of 77? obtained 
here, effectively scaling each of the two sensor data curves 
shown here (as each data set relied on a single pattern 
of heating and thermometer sampling, it is likely that 
each data set would have a systematic error character- 
ized by a single scale factor). Without a more reliable 
model for calculating the thermal integral from the lim- 
ited number of thermometers on the sensor, the data pre- 
sented use the simple temperature interpolation, with an 
estimated systematic uncertainty of order 50% and 20% 
for the TN2mm and EM4mm sensors, respectively. As 
such, the data for both sensors are consistent with the 
predictions for the radiometric effect. The uncertainties 
are however too large to test the radiometric model at 
a level approaching the statistical resolution of the mea- 
surements, or to test definitively the modifications to the 
infinite plate model predicted by the simulations for the 
two sensors. 

In addition to confirming the model of the radiometric 
effect, these data confirm, and extend - to a wider range 
of pressure, temperature, and to a second prototype sen- 
sor - our preliminary observations with the TN2mm sen- 
sor [1^, [l3|. We note that a thermal gradient-induced 
torque that increases with pressure is also observed in 
a LISA-related study at the University of Washington 
(UW) [2l[, for measurements at two pressures and two 
temperatures. This study, which employed a much sim- 
pler and more open geometry than the LISA GRS, ob- 
served forces of the same order of magnitude as those ob- 
served here. However the authors report an unexpected 
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FIG. 10; Normalized measured slopes 4ro-9p^ of the fit- 
ting lines in Fig[7] as a function of the average tempera- 
tures To (blue dots: prototype TN2mm; red dots: prototype 
EM4mm) . The results represent measurements of the radio- 
metric correction factor jr, which is unity in the infinite plate 
model. Horizontal dashed lines: values predicted by the sim- 
ulation for a linear thermal gradient (blue and red refer to 
TN2mm and EM4mm, respectively) from App. |X] We note 
that the error bars represent only the statistical uncertainties 
in the torque measurement. The systematic uncertainty in 
the estimate of the relevant thermal integral is much larger, 
of order 20% and 50% for the EM4mm and TN2mm proto- 
types, and serves as an uncertainty in the overall scale factor 
for each of the two datasets. 



frequency dependence, which is not observed at all in 
the measurements presented here. In addition, the UW 
study reports an increase in the pressure dependent ef- 
fect with increasing temperature, which is contrary to 
the expected radiometric dependence observed here. 
Given the differences in the experimental geometry and 
the limited number of presented measurements, it is not 
possible at this time to further discuss these apparent 
discrepancies. 



B. The zero pressure limit: thermal radiation 
pressure and the effect of outgassing 
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FIG. 11: Measured intercepts of the linear fittings of the 
curves reported in Fig[7l Top panel: data of the TN2mm 
prototype, represented with blue dots. Bottom panel: data 
of the EM4mm, represented with red dots. Error bars indi- 
cate only the statistical uncertainty in the torque measure- 
ments. The lines of the radiation pressure contribution to the 
intercepts predicted by different models are shown for com- 
parison in both panels. The brown continuous lines show the 
prediction of the infinite plates approximation (ffTg as in 
Ean lll|l . which is the same for both sensors. The other mod- 
els are derived from the numerical simulations discussed in 
App. |B]and using the temperature patterns estimated with 
the procedure described in App. [C] Different values have 
been taken into account for the surface absorptivity a, and 
we indicate with d — 100% the case of diffuse light scattering 
and with d — 0% the specular light scattering case. Orange 
line: a = 100%, d = 0%. Black: a = 10%, d = 100%. Ma- 
genta: a = 10%,d = 0%. Green: a = 5%, d = 100%. Cian: 
a = 10%, d = 0%. 



In Fig llll we plot the zero-pressure values of the nor- 
malized torques, N*{P — 0), extrapolated to P = from 
the data in Fig[71 as a function of the average tempera- 
ture for the two sensors. Shown for comparison are dif- 
ferent estimates of the radiation pressure contribution. 
These estimates, evaluated for the experimental temper- 
ature distributions as estimated in Appendix C, include 
the infinite plate model of Eqn[Tl] (brown), and the re- 
sults of the radiation pressure simulation presented in 
Sec. |B]) for 5% absorption and both diffuse (green) and 
specular (cian) reflection. As discussed in Sec. |B1 for 



the geometries of these GRS prototypes, the expected 
thermal radiation effect is substantially suppressed, rel- 
ative to the infinite plate prediction, a suppression that 
is enhanced with increasing reflectivity. Reflectivities of 
95 % (a = 0.05) or higher are in fact expected for gold 
coating at reflectivity thermal radiation wavelengths [27| . 
While the data for the TN2mm are consistent with the 
presence of the radiation pressure effect, the data for the 
EM4mm prototype dip slightly below even the 95% re- 
flection prediction for the radiation pressure effect at the 
lowest temperatures. This discrepancy will be addressed 
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FIG. 12: Two measurements of the thermal induced normal- 
ized torque for the EM4mm sensor as a function of pressure, 
at roughly the same average temperature (blue 298.9 K, red 
299.7 K) After the first measurement (blue), the sensor tem- 
perature was progressively increased to 314 K (highest tem- 
perature in FigO while measuring at the higher temperatures 
of the study. The second (red) data set was measured subse- 
quently after cooling to 299.7 K. The slopes agree within the 
measurement uncertainty after correction for the small tem- 
perature difference. The measured intercept, on the other 
hand, is reduced by more than a factor 2 in the second mea- 
surement, indicating a sharp reduction of the out gassing effect 
during the "bakeout" at 314 K. 



shortly. 

For both sensors, we observe an increase in the zero- 
pressure torque intercepts with temperature, as ex- 
pected both for the radiation pressure and outgassing 
effects. However, for both sensors, the observed increase 
with temperature is clearly faster than the depen- 
dence expected for radiation pressure, regardless of the 
surface reflection properties or geometry, and suggests 
the exponential-like dependence expected for outgassing. 
Further evidence for the existence of outgassing-related 
phenomena is found in the hysteretic behavior of the in- 
tercept following a thermal cycle (see Fig fT2|) . This shows 
that even a mild "bakeout" at 314 K suppresses the inter- 
cept near 299 K by more than a factor two, a hysteresis 
that must reasonably be associated with the degassing of 
the surface at warmer temperatures. 

While outgassing is thus the main effect contributing 
to the zero pressure data, we note that the its magnitude 
near 293 K is, for the TN2mm prototype, only roughly 
equal to that calculated for the infinite plate model for 
the radiation pressure effect and, for the EM4mm sensor 
designed for LISA Pathfinder, considerably smaller. As 
such, the outgassing effect as measured here does not 
pose a significant increase to the thermal gradient noise 
budget already considered for LISA. 

Even if radiation pressure plays a negligible role in the 



zero-pressure torque, no simple model can explain the 
observation that, for the EM4mm prototype, the extrap- 
olated zero-pressure intercept of the normalized torque 
becomes negative at the lowest temperatures studied. 
This small - roughly 15% of the magnitude of the in- 
finite plate prediction for the radiation pressure effect - 
but statistically significant violation of the expected sign 
for a thermal gradient induced torque is visible in the 
283 K data point in Fig[TlJ There are several possible 
explanations for this observation: 

• Pressure offset If the ion gauge pressure measure- 
ments were higher, by a constant offset Pqff, than 
the actual pressure inside the sensor, then the es- 
timated zero-pressure torque intercepts would be 
artificially lowered. This hypothesis is at least con- 
sistent with the fact that, in spite of the fitted neg- 
ative zero-pressure torque intercepts, we do not ob- 
serve negative torques even at the lowest measured 
pressure values. Pqff would have to be roughly 
2 X 10~^ Pa to explain the observed negative torque 
offset. Given the quoted accuracy of the Pfeiffer 
gauge, this offset should not exceed 1.5 x 10~^ Pa 
(15% of the 10^^ Pa full scale used for the low 
pressure data). Such a pressure measurement offset 
could thus explain most, but not all, of the observed 
negative torque intercept. Another, slightly less 
stringent, upper limit on Pqff comes from the low- 
est pressure measured in the study, 2.3 x 10~^ Pa. 

We also note that an offset arises in the outgassing 
from inside the sensor, which causes the true pres- 
sure inside the sensor to exceed that measured in 
the vacuum chamber. This thus has the wrong 
sign to explain the measured negative zero-pressure 
torques. Additionally, this effect is estimated to be 
quite small; given the effective impedance of the 
holes connecting the sensor to the rest of the vac- 
uum chamber, the pressure inside the sensor would 
be only 10^^ Pa higher than outside, even if as 
much as 10% of the total measured apparatus out- 
gassing (roughly 3 x 10~^ mbar 1 / s) originates 
inside the GRS. 

• Inhomogeneous outgassing Outgassing, associated 
with contamination and "virtual leaks," such as 
the joints between different pieces, is likely to be 
highly dependent on position inside the GRS. It 
is possible that, at least for certain temperatures, 
the GRS outgassing could be dominated by one or 
several localized points that are in positions, likely 
near the corners of the sensor, that are warmed in 
phase with the driven heaters during the experi- 
ment, but nonetheless contribute a torque of the 
opposing sign. 

At present we can not discriminate between these or 
other explanations for the negative zero-pressure torque 
observed at 283 K for the EM4mm prototype. 
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V. CONCLUSIONS 

This paper represents the first extensive study aimed 
at characterizing thermal-gradient induced forces in the 
specific environment relevant to LISA free-fall, realisti- 
cally capturing the key environmental parameters that 
the LISA test masses will experience inside of the gravi- 
tational reference sensors in orbit. This includes not only 
temperature and residual gas pressure, but also geometry, 
materials, construction technique, and handling history. 
The conclusions of the study, based on torsion pendulum 
measurements made with two representative GRS proto- 
types and numerical simulations of the radiometric and 
radiation pressure effects, is that thermal gradient forces 
will likely not threaten the LISA free-fall goals. 

For the radiometric effect, the calculated thermal gra- 
dient transfer function is roughly 25% larger than that 
predicted by the simple infinite plate formula. This in- 
crease - yielding |^ « 23pN/K or « 1.25 - is due 
largely to the shear forces exerted by molecules strik- 
ing the Y and Z faces of the TM. In the experiments 
performed, the radiometric effect can be isolated from 
the pressure dependence of the measured torques in the 
measurements present here. The expected functional 
dependence, as well as a frequency-independent linear re- 
sponse, is observed. The amplitude is quantitatively con- 
sistent with that predicted by the simulations, for both 
sensor prototypes, but the systematic uncertainty in the 
estimation of the relevant temperature profile - estimated 
to be roughly 20% for the prototype more closely repre- 
senting the current LISA design (EM4mm) - is not suf- 
ficient to distinguish the corrections to the infinite plate 
model foreseen in the simulations. 

/^From the radiation pressure simulations, we observe 
how larger gaps allow a dilution of the force from bounc- 
ing thermal photons, given the highly reflective GRS en- 
visioned for LISA. Assuming 95 % reflectivity and spec- 
ular reflection, the resulting radiation pressure trans- 
fer function is reduced to roughly a third of the in- 
finite plate model prediction (k^p w 0.32), yielding 
IStIb^p ^ 9pN/K. Radiation pressure, along with out- 
gassing and any other pressure-independent effect, are 
studied with the torque data in the limit of zero pres- 
sure. These data here are dominated by the outgassing 
effect, clearly distinguished from the radiation effect by 
the sharp temperature dependence and the observed re- 
duction in the effect following a mild "bakeout." How- 
ever, the effect is not large enough to threaten the LISA 
noise budgets; for both prototypes, the torque at room 
temperature is comparable to (the older 2 mm-gap pro- 
totype sensor) or smaller than (LISA 4 mm design) that 
foreseen in the infinite parallel plate model for radiation 
pressure. Additionally, it is likely that the outgassing 
effect, already tolerable for LISA in these sensors which 
have never had a serious bake-out, would very likely be 
further reduced with the 180 K bakeout foreseen before 
flight. 

The measured torques are thus consistent with the pre- 



dicted radiometric effect and reveal no large zero-pressure 
effect, from outgassing or anything else, that exceeds the 
levels already included in the previously-assumed radi- 
ation pressure estimates. As such, these experiments 
increase confidence in the LISA noise model and sug- 
gest that the overall thermal gradient transfer function 
is no larger than that indicated in the introduction, 
2^ ~ lOOpN/K. One caveat, however, is that these 
experiments are based on torque measurements, rather 
than the x force relevant to LISA. Outgassing is not likely 
to be homogeneously distributed within the sensor, and 
the torque measurements are sensitive to most, but not 
all outgassing locations, particularly those acting cen- 
trally on the TM face. To address this concern, we are 
currently developing a new torsion pendulum, with the 
TM displaced from the torsion axis, that will be directly 
sensitive to the thermal gradient-induced x forces [32| . In 
addition to a direct measurement of the quantity of inter- 
est for LISA free-fall, the pendulum geometry, sensitive 
to the translational temperature difference AT^ will al- 
low a much simpler analysis of the relevant temperature 
profile. 



APPENDIX A: RADIOMETER EFFECT: 
NUMERICAL SIMULATIONS DETAILS 

The basic idea is to simulate a gas of non-interacting 
particles bouncing in the gaps between the TM and the 
EH, in the presence of a thermal gradient along the EH. 
The ballistic motion of a single particle was calculated, 
with the impulse and moment of the impulse recorded 
at every collision of the molecule with the TM. The to- 
tal TM force and torque exerted by the particle are cal- 
culated by multiplying by the number of particles that 
would fill the gaps at a given mean temperature and pres- 
sure. 

The sensor geometry considered is that of a cubic 46 
mm TM placed in the center of a rectangular box, both 
of which are fixed during the simulation. There is no 
hole in the sensor that would allow exchange of parti- 
cles with the surroundings. As edge effects are expected 
to be more important when the gap between the TM 
and the EH increases, several configurations with differ- 
ent gaps has been analyzed, including the design under 
construction for LISA Pathfinder and currently proposed 
for LISA (this geometry is represented in Fig[T3] by the 
point at 3.2 mm, which is the average of the gaps along 
the Y and Z axes, 2.9 mm and 3.5 mm, respectively). 
The gap sizes considered range from 1 mm to 8 mm. 

The calculations for the translational force for LISA 
are performed with a linear gradient across the position 
sensor along the x axis, with the inner X faces Xj^, X- 
each at uniform temperature and the inner surfaces of the 
Y and Z faces having a linear temperature profile. The 
TM is assumed to be isothermal at the average tempera- 
ture of the position sensor. We note that the calculation, 
as it has been set up here, must be performed with a 



temperature distribution defined in every point in the 
sensor, and any modification of the temperature distri- 
bution, even locally, must be studied as an entirely new 
simulation. 

At every bounce onto a surface, the particle is assumed 
to thermalize with the surface, with reemission that has 
no memory of what happened before the bounce. When 
the particle leaves a surface, we extract the velocity direc- 
tion from a Knudsen distribution [3l| , with flat probabil- 
ity both in the azimuthal angle (j> and in cos^(0), where 
is the angle with respect to the surface normal vec- 
tor. The kinetic energy is extracted from a distribution 
of velocities that would produce Maxwell-Boltzmann dis- 
tribution in the space. 

Each time a particle hits a TM face, it completely 
transfers its momentum to the TM. Upon reemission, it 
transfers the recoil momentum to the TM and then flies 
at constant velocity until it reaches another surface. Mo- 
menta and angular momenta transferred to the TM are 
summed bounce by bounce, as well as the flying time be- 
tween bounces. At the end, the total linear and angular 
momenta transferred to the TM are divided by the total 
elapsed time, giving the average force and torque acting 
on the TM. This calculation has been performed for each 
momenta component transferred to each TM face. For 
each geometrical configuration, roughly 10 runs each con- 
sisting of roughly 5 x 10^ bounces were performed. The 
resulting force and torque values were averaged and the 
standard deviation calculated to determine the calcula- 
tion uncertainties. 

The results of the calculation are summarized in 
FigHniThe GRS design under construction for LISA 
Pathfinder and currently proposed for LISA, is repre- 
sented in Fig[T3]by the point at 3.2 mm. In Sec. Ill At hese 
results are illustrated and discussed. 



APPENDIX B: THERMAL RADIATION 
PRESSURE: NUMERICAL SIMULATIONS 

The effective role of thermal radiation pressure has 
been evaluated in the simplified GRS geometry, with a 
cubic TM contained and centered inside a rectangular 
box, by numerical simulations. The eventual modifica- 
tions to the radiation pressure "transfer function" can 
be calculated and expressed as the correction factors khp 
and jjip introduced in Sec. Ill Bl for comparison to the 
simple model discussed there. 

For the simulations, the total forces and torques on the 
TM associated with thermal radiation have been repre- 
sented as an integral of multiplied by an independent 
vectorial force or torque function, ^{r) or lt{r), rep- 
resenting the contribution (per K^) of thermal photons 
originating from the sensor surface element ds at position 
r of the EH surfaces S : 

^RP = E / 7(r)r4(r)ds 

(EH.TM faces) 
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FIG. 13: Results of numerical simulations for the radiometer 
effect, as a function of the gap (in mm) between test mass 
and electrodes of the position sensor. The top panel shows 
the ratio of the simulation torque to that calculated from 
EqnO for the x component of the force acting on the TM X 
faces (red dots), for the x force acting on the Y and Z TM 
faces (blue dots) and for the total x force acting on the TM 
(sum of previous values, shown as green dots). The bottom 
panel shows the ratio of the simulated torque, around the z 
(or y axis, to the infinite plate result (Eqn[5ll , for the Y (or Z) 
face. The average of the magnitudes for the 4 X — Y faces has 
been reported. While all other calculations are made for equal 
gaps in all directions, the point at 3.2 mm is calculated with 
the real gaps of the configuration baselined both for LISA 
Pathfinder and LISA and shown in Fig[l] 



(EH, TM faces) 

In the infinite plate limit (see Eqn[5]), fx — for sen- 
sor surface elements directly opposite the TM faces on 
the X+ and X- faces (the yellow and orange zones at 
left in Figin]), respectively, and zero elsewhere in the sen- 
sor. For the torque, 6(s), where the effective 
armlength b(s) = ±y for the orange and yellow zones 
at right in FiglS] (and, analogously, ±x for the corre- 
sponding zones on the Y faces). The functions ~f(r) and 
Tf (r) depend only on the sensor geometry and reflection 
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properties, and so can be evaluated point-by-point and 
independently of the sensor temperature distribution. 

In the simulations, a chosen number of photons is emit- 
ted from a given position of the TM or of the EH surfaces. 
Photon angle is determined statistically, again assuming 
a uniform probability density in cj) and cos^ (0). The sur- 
face absorption properties were assumed to be indepen- 
dent of wavelength, which allowed for a simple normaliza- 
tion for the number of photons (or effective time of the 
simulation) and the resulting forces, using the Stefan- 
Boltzmann law for the radiated power from a surface 
element of given absorption coefficient a, (acrT*) , and 
the energy - momentum relation for photons, E — pc. 
Each emitted photon is propagated ballistically from the 
emitting surface, and then either absorbed or reflected 
statistically, according to the surface absorption proper- 
ties imposed. In the case of reflection, specular or diffuse 
(statistical angular reemission) was decided statistically 
with the imposed percentage of diffuse scattering, d. At 
each collision with the TM, the momentum and moment 
of the momentum transfer are recorded and then summed 
at the end of the simulation to determine the total force 
and torque on the TM. 

For the cases studied here, emission from the TM was 
not calculated, as, in the isothermal TM-approximation 
reasonable for LISA and our experiments, the TM radia- 
tion produces no net recoil. Additionally, for simplicity, 
the absorption properties were assumed to be uniform 
and equal for the TM and sensor surfaces. The points 
of emission for which the simulation was performed were 
uniformly distributed on the faces of the sensor, using a 
rectangular grid with 0.5 mm spacing. Each simulation, 
typically involving 200 photons, was repeated, typically 5 
times, for each point inside the sensor, in order to obtain 
the uncertainties for ~f{r) and 7f (r). 

Simulations have been performed for several different 
sensor geometries, including the two studied experimen- 
tally in this paper. Additionally, simulations have been 
performed for several sets of reflection properties: (i) 
absorptivity a = 100% (ii) absorptivity — a — 10% 
and pure diffusive reflection (d = 100%) (iii) absorp- 
tivity a = 10% and pure specular reflection {d — 0%) 
(iv) absorptivity — a = 5% and pure diffusive reflection 
(d — 100%) (v) absorptivity a = 5% and pure specular 
reflection (d = 0%). 

Finally, in order to obtain the relevant radiometric ef- 
fect correction factors for the x-force and the x-torque, 
krp and 7flp, the values of 7(r) and Tf (r) obtained from 
the simulation must be multiphed by (T(r))^, with the 
sensor temperature distribution T(r) chosen to represent 
a specific experimental thermal perturbation. Results for 
cases relevant to thermal gradient forces and torques are 
discussed below. 



1. Radiation pressure simulation results: induced 
force 

For the study of the radiation pressure force along x, 
we plot first the relevant force function /x(r), normalized 
by the infinite plate value {2a /3c), as a function of the 
emission position r on the sensor X and Y faces (the Z 
faces play a role very analogous to that of the Y faces for 
the X force). These are shown for the LISA Pathfinder 
sensor design in Fig[T3] {X faces - left column, Y faces - 
right column), together with infinite plates model case. 
Then, to calculate the relevant factor kad, integration 
of fx is performed for a temperature distribution with a 
gradient in the x direction, with each X face at a uniform 
temperature and a linear temperature gradient across the 

Y and Z faces. The resulting values for the correction 
factor Kflp are summarized for different sensor geome- 
tries and absorptivity characteristics in Table [Til The ta- 
ble also gives a breakdown of the contributions to khp, 
including not only the contribution expected from the X 
faces directly opposite the TM {kx{in)): but also those 
not foreseen by the infinite plate model, from the bor- 
ders of the X sensor faces {kx(out)) and from the other 
sensor faces (kyz)- 

As visible in FiglT4j2"'' row, left), where we con- 
sider full absorption of emitted photons (a = 100%), 
for emission from near the center of the sensor X faces, 
fx/ (2ct/3c) ~ 1, as expected for the infinite plate model 
(1** row left). Approaching and crossing the edge of the 
TM, fx/ (2cr/3c) decreases smoothly from one towards 
zero, as a continuously smaller fraction of the emitted 
photons strike the TM surface. For runs with the small- 
est gaps (not plotted here), the transition near the edge 
of the TM becomes sharper, and the infinite plate model 
becomes more accurate. Additionally, there is a signif- 
icant contribution of emission from the Y faces to fx 
Fig fTW right column). The zones near the edges adjacent 
to the A"+ and X^ faces contribute photons that are pre- 
dominantly absorbed on the TM near the center of the 

Y face, and as such exert a shear force along x. There is 
an analogous contribution from the Z face. When inte- 
grated with the temperature distribution, this shear force 
contribution is largely responsible for the increase in the 
radiation pressure effect {kup — 1-17, Table HIl). We note 
that this increase in the case of full absorption is simi- 
lar to that seen for the radiometric effect, with the shear 
forces along the Y and Z faces a leading factor in both 
cases. Also in Table HIl we see that with 200 micron gaps, 
the simulation confirms the infinite plate model, krp = 1 
to better than 1%, nearly entirely from the contribution 
of the X faces opposite the sensor Kx{in)- 

With higher reflectivity, i.e. low a (3'"'^ to 6*'' row, left), 
the effect of smoothing the X face contribution is more 
pronounced, as even photons emitted near the center of 
the X face can migrate off the TM X face and deposit 
their momentum on other faces, even on the opposite TM 
X face. The smoothing is more pronounced for specular 
scattering than diffuse, because in diffuse scattering the 
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FIG. 14: Intensity plots of the x component of the force func- 
tion fx{r), the force per K* per unit area of emitting sensor 
surface if*, normalized by {2a /3c), the characteristic value of 
the infinite model. Results are shown for emitting points on 
the X (left column) and Y surfaces (right column), with the 
outline of the TM section a dashed-line 46 mm square in each 
case. The 1"* row shows the data for the infinite plates model. 
Rows 2 to 6 show the different cases of refiection properties: 
full absorption, a — 100%, absorption coefficient a = 10% or 
a = 5%, pure diffusive reffection (d = 100%) or purely spec- 
ular refiection (d — 0%). Legends for the intensity values are 
shown in the colors bars on the right of each plot. 



GRS Geometry 
(TM , gap) 


(absortivity, 
diffusion) 


KflP Kx(IN) I^X{OUT) I^YZ 


(46mm^, 4mm) 


(100%,0%) 
(10%, 100%) 

(10%, 0%) 
(5%, 100%) 
(5%, 0%) 


1.17 0.95 0.05 0.17 
1.01 0.69 0.05 0.28 
0.63 0.50 0.02 0.10 
0.75 0.49 0.03 0.23 
0.32 0.26 0.01 0.05 


(40mm^, 2mm) 


(100%,0%) 
(10%, 100%) 
(10%, 0%) 
(5%, 100%) 
(5%, 0%) 


1.12 0.96 0.04 0.12 
1.01 0.78 0.03 0.21 
0.75 0.64 0.02 0.10 
0.85 0.62 0.03 0.20 
0.45 0.39 0.01 0.05 


(46mm^, 0.2mm) 


(100%,0%) 
(10%, 100%) 
(10%, 0%) 


0.999 0.986 0.002 0.01 
0.966 0.937 0.002 0.027 
0.946 0.906 0.002 0.038 


(40mm^, 0.2mm) 


(100%,0%) 
(10%, 100%) 
(10%, 0%) 


100 0.985 0.002 0.013 
0.963 0.929 0.002 0.031 
0.942 0.896 0.002 0.043 



TABLE II: The overall correction factor hrp estimated in the 
case of a linear temperature profile across the sensor along x 
axis, and the X-|- and X- faces at uniform temperature;^ = 
{t^xi^ + Ka;„„t -I- and t^x^^t a-^d K{j/-fz) represent- 

ing respectively the contributions from the X faces facing 
the TM, the X faces not facing the TM and the lateral Y 
and Z faces. The simulations were performed for different 
GRS geometries (i.e. different gap); the one indicated with 
(46mm^, 4mm) is the GRS geometry under construction for 
LISA Pathfinder and currently baselined for LISA. Statistical 
errors in each estimate of nnp axe in the order of 2 4%. 



effective distance traveled by the photons increases only 
as the square root of the number of bounces, in a random 
walk, rather than hnearly as for specular reflection. The 
extreme case of high reflectivity and specular reflection 
(a = 5% and d = 0%, bottom left) represents a pro- 
gressive homogenization of the radiation pressure, as a 
photon generated from any point in the sensor winds up 
imparting momentum on different parts of the TM, effec- 
tively washing out also the shear effect (see FigfHl bot- 
tom right). As a value for adsorption around 5 -=-10% and 
specular reflection is the likely case for the gold coated 
surfaces envisioned for the LISA GRS, a reasonable esti- 
mate of the radiation pressure correction for LISA ranges 
between krp = (0.32 ±0.01) and krp = (0.63 ± 0.02) 
(the errors given are simply the statistical uncertainties 
of the simulation). 



2. Radiation pressure simulation results: induced 
torques 

In Fig llSI we show the z components ^^(r) of the 
torques per unit area per K** as a function of the position 
inside the EH and as a function of the GRS surface prop- 
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FIG. 15: Intensity plot of the torque per unit area (r) about 
the z axis, normalized here by the factor {2a /3c), for emitting 
points on the X surface (the y face of the EH shows the exact 
same profile), the 1** plot (top left) shows the infinite plates 
model case. The various cases of full absorption, a — 100%, 
absorption coefficient a — 10% or a = 5%, pure diffusive 
reflection (d — 100%) or purely specular reflection (d — 0%) 
are shown. 



erties. We still use the normalization factor (2cr/3c): we 
thus remark that nz{r)/ {2a/3c) ~ 6(r) for the infinite 
plate approximation (shown in Fig|15l top left), which 
coincides in the particular case of the X face of the EH 
under consideration with the y coordinate, and with the 
X coordinate on the Y face. The different cases of surface 
reflectivity are here shown (left to right, top to bottom). 

As is visible, a reduction similar to that observed for 
the force is seen for the thermal radiation torque. How- 
ever, in this case, the effect is much more relevant: the 
overall attenuation of is much more effective, because 
the attenuation of /z(r) is progressively more effective 
approaching the edges of the TM. It thus gets reduced ex- 
actly where the arm-length b{r) is bigger and the torque 
per unit area should contribute most to the overall ef- 
fect. This is clearly visible in the data in Table [1111 where 
we report the estimates of the radiation pressure induce 
torque correction factors 7i?p, for an ideal linear temper- 
ature profile on one of the faces of the IS (with zero tem- 
perature modulation on the other faces) . We estimate the 
total torque to be reduced down to « 73% compared to 
the theoretical estimate in the case of 100% absortivity. 
In the case of absortivity of the IS surfaces of order 10% 



or lower, the reduction is then much more relevant, and 
we estimate the total torque to be suppressed down to a 
value between 15% for EM4mm and 24% for TN2mm of 
the theoretical value, in the realistic case of reflectivity 
95% (a = 5%) and pure specular reflection {d = 0%)(see 
Table Iml). 



GRS Geometry 


(absorptivity. 


7tot 


(TM , gap) 


diffusion) 






(100%,0%) 


0.73 


(46mm^, 4mm) 


(10%, 100%) 


0.30 




(10%, 0%) 


0.27 




(5%, 100%) 


0.17 




(5%, 0%) 


0.15 




(100%,0%) 


0.80 


(40mm^, 2mm) 


(10%, 100%) 


0.42 




(10%, 0%) 


0.39 




(5%, 100%) 


0.27 




(5%, 0%) 


0.24 



TABLE III: Correction factors 7_rp for the thermal radiation 
pressure induced torques. The GRS geometries here reported 
(46mm^, 4mm and 40mm 2mm) correspond to the two GRS 
prototypes studied with the torsion pendulum (EM4mm and 
TN2mm respectively). Statistical errors are in the order of 
2 ~ 4%. 



APPENDIX C: ESTIMATE OF TEMPERATURE 
PATTERNS AND THERMAL INTEGRALS 



As mentioned in Sec. IIII CI for each measurements 
cycle we evaluate the component at frequency lu of the 
"thermal torque" (/^ T(s)6(s)ds) with the equation 

/ Tis)b{s)ds) = E / i±)T{y,z,u;)ydydz + 

J S / uj Xj^.X- 

+ J2 [ iT)T{x,z,uj)xdydiCl) 

The components T(si, S2,i-^) of the temperature at 
position(si, S2), where si and S2 are two coordinates 
mapping the 4 X-Y sensor surfaces, have been esti- 
mated by a cubic spatial interpolation of the components 
T{si, S2,uj)i of the readings Ti{t) of the thermometers 
on the external electrode housing surface at the posi- 
tions (si, 52)1, which were measured during assembly. In 
performing the interpolation, several important approx- 
imations were made and are addressed in the next two 
paragraphs. 

Figure [5] shows a cartoon of the two GRS proto- 
types with the locations of the thermometers (small blue 
squares). For the EM4mm prototype, six of the eight 
installed thermometers were used for the interpolation, 
namely Tl, T2, T3, TA, T7, and T8. For TN2mm, four 



of the five thermometers were used for the interpolation 
(r2, TS, T5 and T6). In this case, the reading of T5 was 
assigned to the temperature at the center of the heater 
HIE, in order to constrain the local temperature maxi- 
mum to coincide with the heater location. Additionally, 
in order to construct the temperature distribution on the 
TN2mm X+ face, which has only a single thermometer, 
in contrast with the three X— face thermometers, we in- 
troduced an artificial thermometer reading, Tfake near 
the heater H2W. The modulation of Tfake was assigned 
to be the value of the interpolated component at the cen- 
ter of the electrode below H2E on the opposing X— face, 
scaled by the ratio in order to account for observed 
asymmetric heating for the two X faces. 

Figure [TBI shows an example of a reconstructed temper- 
ature profile T(si, S2,u}) across the entire sensor TN2mm. 
Several assumptions were made in making this estima- 
tion, for both sensors. First, the temperature on the 
inner surface of the sensor is assumed to be equal to that 
on the corresponding outer surface position, as calculated 
by the interpolated the outer surface thermometer read- 
ings. Additionally, the temperature profile along the z 
axis (that parallel to the torsion fiber in our experiments) 
was collapsed to a single value, that of the temperature 
readings interpolated in the horizontal dimension. The 
measured temperatures are thus assumed to represent the 
average temperatures along the sensor z axis at a given 
value of the horizontal coordinate x (or y, depending on 
the face). These assumptions clearly result in some sys- 
tematic error in the estimation of the thermal integral, 
which is discussed next. 

The temperature read by a thermometer on the GRS 
external surface is a good approximation of that mea- 
sured on corresponding the inner GRS surface position, 
as long as the thermometer is not too close to an oper- 
ating heater or to relevant localized thermal interfaces, 
such as the supporting legs at the four lower corners of 
the sensor. All the thermometers used in the interpo- 
lation are indeed far from these corner supporting legs, 
and, additionally, they are all placed at least 9 mm from 
the border of the heaters, a distance equal to the sen- 
sor wall thickness, including the electrode. A simplified 
FEM model has confirmed that the inside-outside tem- 
perature difference can indeed be considered negligible in 
these conditions. 

As stated above, the thermometer readings used for 
the interpolation are assumed to be the average temper- 
ature value along z. For the Y faces, where there are 
no heaters, the main effect of the temperature variation 
along z is expected to be an altitude-dependent offset of 
the whole profile, with a negligible impact on the ther- 
mal integral T{s)b{s)ds, where a constant temperature 
offset integrates to zero. For the X faces, however, the 
operating heaters produce a more complicated and more 
rapidly varying temperature pattern. Moreover, the in- 
ner surface temperature modulation in correspondence 
with the operating heater is expected to be relatively 
uniform across the heater, but higher than that mea- 




) OS 12 

position s [m] 

FIG. 16: An example of the temperature profile across the 
entire sensor TN2mm. Blue line: derived from cubic inter- 
polation of the thermometer readings. Red line: based on 
a constant temperature across the electrodes with operating 
heaters and a linear interpolation of the thermometer read- 
ings across the rest of the sensor. A single temperature is 
assigned for all values of z at each given value the horizontal 
coordinate s (x or y, depending on the face). 



sured by the closest thermometer. A rough estimation 
of the consequent correction to the temperature pattern 
on the X faces has been obtained by assuming that the 
heat flux in the proximity of a cylindrical heater is ba- 
sically radial and uniform, an approximation supported 
by a simplified FEM calculation. In this approximation, 
we estimate that our simple interpolation technique un- 
derestimates the thermal integral by roughly 30% for the 
UTN2mm and 20% for the EM4mm. However, given the 
approximate nature of the analysis and the experimen- 
tal difficulty of verifying this model with a tight array 
of thermometers on the inner sensor surfaces, we choose 
here not to correct the data for this likely systematic 
error. Instead, we present the data obtained from the 
relatively straightforward interpolation and discuss the 
role of this systematic error in the conclusion of our data 
analysis. 

In order to address the dependence of the thermal in- 
tegral calculation on the method of spatial interpolation, 
we compare the results of two different methods, both of 
which reasonably reproduce the expected relatively uni- 
form temperature profile of the "hot" electrodes coin- 
ciding with the locations of active heaters (see Fig fTH)) . 
The first calculation assumes a perfectly constant tem- 
perature across the "hot" electrodes and a linear inter- 
polation for the rest of the temperature profile. These 
results are compared with that obtained with the cubic 
polynomial interpolation. Thanks to the higher number 
of thermometers, for the EM4mm prototype the differ- 
ence was negligible, while for the TN2mm model the first 
method gave a value 20% larger than the second. 

Considering the approximate analysis of the different 
sources of uncertainty in the thermal integral, we conser- 
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vatively assume the overall uncertainty to be the linear 
sum of the above estimated errors. This yields an un- 
certainty of roughly 50% for the TN2mm prototype and 
roughly 20% for EM4mm. Given that the same heat- 
ing pattern is repeated in all measurements presented for 
each of the two sensors, this uncertainty in the thermal 
integral estimation can be thought of as a single scale 
factor for the each of the two sensor data sets, rather 
than an uncertainty acting incoherently on different data 
points in the same sensor. Unless otherwise specified, 
the error bars of the torque data normalized by means 
of the "thermal integral", namely N* represents just the 
statistical error of the torque measurement N. The role 



of this scale factor will be discussed with the conclusion 
of our data analysis in Sec. IIVI 
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